Chaos synchronization in a hyperbolic dynamical system with long-range interactions 
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We show that the threshold of complete synchronization in a lattice of coupled non-smooth chaotic 
maps is determined by linear stability along the directions transversal to the synchronization sub- 
space. As a result, the numerically determined synchronization threshold agree with the analytical 
results previously obtained [C. Anteneodo et al, Phys. Rev. E. 68, 045202(R) (2003)] for this class 
of systems. We present both careful numerical experiments and a rigorous mathematical explanation 
confirming this fact, allowing for a generalization involving hyperbolic coupled map lattices. 



The possibility of synchronizing chaotic dynamics has 
been harnessed in a large number of systems ofphysi- 
cal interest [![, like coupled Josephson junctions 0] and 
lasers [||. Although there have been identified different 
types of chaos synchronization, we shall concentrate on 
the so-called amplitude or complete synchronization, for 
which all dynamical variables undergo the same time evo- 
lution Q • The essential dynamics involved in the process 
of chaos synchronization lies on the low-dimensionality of 
the subspace (in the phase space of the system) in which 
synchronized motion sets in. 

For example, if we consider a lattice of N coupled os- 
cillators, each of them represented by a vector field of D 
dimensions, where typically D <§; N, the synchronized 
state belongs to a D-dimensional subspace of the ND- 
dimensional phase space. In order for this synchronized 
state to exist the coupling among oscillators takes on a 
suitable form [5| . Whether or not this synchronized state 
is stable, however, is a more difficult question, since it in- 
volves the analysis of infinitesimal displacements from the 
synchronized state along all (N — 1)D directions transver- 
sal to the synchronization subspace Q . The stability con- 
dition of the synchronized orbit with respect to transver- 
sal perturbations can be obtained from the negativeness 
of the largest transversal Lyapunov exponent. 

In this letter we consider a coupled chaotic map lat- 
tice (CML) in which the coupling prescription is non- 
local, for it takes into account the distance between them 
along the lattice. Such non-local couplings appear in 
many problems of physical 0] and biological interest § . 
We suppose that the coupling strength decreases with 
the lattice distance as a power-law, which characteris- 
tic exponent can take on any positive value 0]. The 
loss of transversal stability of the synchronized state, as 
the coupling parameters are varied, was found in such 
power-law couplings, with help of the largest transversal 
Lyapunov exponent, for a number of chaotic maps [l(J. 
In the particular case of maps with constant eigenvalues 
of the Jacobian matrix (piecewise-linear chaotic maps) 
we obtained analytical results for the loss of transversal 
stability of the synchronized state which agree with the 
numerical simulations Such CML's represent hyper- 
bolic dynamical systems, what enables us to use powerful 
mathematical tools like ergodicityand global shadowing 
of numerically generated orbits [121 ] . 



On the other hand, in a recent paper there was ar- 
gued that in the special case of coupled non-smooth maps 
the synchronization transition would not be given by the 
largest transversal exponent, but rather by a different ap- 
proach taking into account finite distances from the syn- 
chronized state [13J . To investigate this apparent contra- 
diction we considered in this letter the transient behavior 
of the non-synchronized orbits for coupled piecewise lin- 
ear maps. Our results show that the analytical results 
of Ref. [ll[ (using linear transversal stability of the syn- 
chronized state) hold for both smooth and non-smooth 
maps, the numerical results being strongly affected by 
many factors as the large transient time and the choice 
of initial conditions. Due to these factors, the time it 
takes to achieve convergence to the synchronized state 
may be extremely large, what may lead to wrong con- 
clusions about the stationary state of the system. Mo- 
tivated by this problem, we investigated the validity of 
the transversal linear stability analysis in a class of hy- 
perbolic CML's, using periodic-orbit theory to unveil the 
role of the unstable orbits embedded in the synchronized 
state HlliEIll- 

The CML we consider in this work can be written in 
the explicit form of a TV-dimensional dynamical system 

x„ +1 = (1 + C)F(x n ) = BF(x„), (1) 

where the components of x„ = (x^P , x^ , ■ ■ ■ , il"') de- 
note the state variable attached to the map located at 
the site i = 1, . . . , N at time n = 0, 1, If the un- 
coupled maps are written as x i— > fix) we can write 

F(x) = [f(x^),f(x^),...J(x^)} T . Moreover, the 
coupling prescription is represented by the matrix C, and 
1 is the identity matrix. 

In the following we consider the generalized Bernoulli 
map f{x) = fix, (mod 1), where x 6 [0,1) and f3 > 1, 
such that the isolated map generates a strongly chaotic 
orbit. When these piecewise-linear maps are coupled ac- 
cording to Eq. (JTJ) , in order to ensure that x„ G [0, 1), 
the elements B must satisfy the following necessary and 
sufficient conditions: By > 0, and < J2f=i Bij < L for 
all i,j = 1,2,... N 0. Moreover, we use a symmetric 
coupling matrix with elements 

Gj=eri- l [ry a Q.-6 iJ )-v8ij], (2) 
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where = mimgzli — j + W| is the minimum lattice 
distance between the sites i and j (with periodic bound- 
ary conditions), rj = 2^f ' r~ a , with N' = (N — l)/2, 
and the coupling strength satisfies < e < 1 due to the 
constraints on _By . The effective range of interactions is 
represented by a > such that the limits a = and 
a — ► oo correspond, respectively, to global (mean field) 
and local (first neighbors) coupling prescriptions. 

A completely synchronized state is the chaotic orbit 
for which x£ = . . . = and which is a solution of 

Eq. ([T]). Since the Jacobian DF = /3B is a circulant 
matrix, its eigenvalues can be analytically obtained as 
AW = /3[(1 - e) + (s/rf)b^-% where 



N> 

& w = y - 



( 2irkm\ 

■ cos I — — — 1 



V n 
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such that the Lyapunov spectrum {Aj} i=1 can be de- 
rived [11]. The stability threshold of the synchronized 
state with respect to infinitesimal transversal displace- 
ments, obtained by imposing A2 = 0, gives two curves 
in the parameter plane (e versus a): (i) s' c (a, N) = 
min{e up (a, AT), 1}; and (ii) e c (a, N) = min{e io (a, N), 1}, 
where we defined 

e up (a,N) = (l+/3- 1 )[l-(6 (Ar %)]" 1 , (4) 
e lo (a,N) = (l-p-^ll-ibW/r,)}- 1 . (5) 

In order to check the validity of these analytical con- 
ditions for the threshold of transversal stability we have 
made careful numerical experiments using the same cri- 
teria as proposed in Ref. [13| (where it has been claimed 
that those conditions would hold only for coupled smooth 
maps). Accordingly, we choose initial conditions Xq uni- 
formly distributed in the interval [0,1) [Hj]. The CML 
is firstly iterated for a transient time of T w = 10 w X N 
times and further iterated by more T = 10 3 x N times. 
As a numerical diagnostic of complete synchronization 
we computed the following quantity 



R = 



E — 

NT 



(0) 



which is essentially a mean deviation from the lattice- 
averaged amplitude. The resulting dynamical state is 
considered as being completely synchronized if R < 10~ 8 . 
In the coupling parameter space we keep a constant and 
sweep through the values of e G [0,1]. The value cor- 
responding to the synchronization threshold, denoted as 
£num, is obtained from bisection as e num — [e s + £,f)/2, 
where e s and 6d are, respectively, the last value corre- 
sponding to a synchronized state and the first value for 
a non-synchronized one. The numerical value of e num is 
turned more accurate from refining the increment mesh 
and repeating the process, until (e s — Ed) < 10~ 3 . 

The results of this numerical procedure, for the case 
(3 = 1.1, are depicted in Figure [TJ where we show the 
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FIG. 1: (color online) Values of the coupling strength at the 
onset of transversal stability loss of the synchronized state, as 
a function of the effective coupling range. We used = 1.1 
and (a) different lattice sizes; (b) different transient times 
T w , for a fixed lattice size; (c) different distributions of initial 
conditions, for N = 257; (d) different initial conditions, for 
N — 33 and a fixed transient time. The solid lines represent 
the analytical results from linear transversal stability of the 
synchronized state. 



value of the coupling strength at the synchronization 
threshold as a function of a. In Fig. HJa) we show how 
the numerically determined critical value increases with 
a for different lattice sizes N, the transient time being 
different for each choice, using w — 5. The solid lines 
correspond to the analytical condition derived in Ref. [Tl| 
(and that depend on the lattice size as well). In fact, as 
the lattice size N increases, the numerical values of e num 
may no longer match the analytically predicted values, if 
a is large enough [l3j . This does not mean, however, that 
the analytical value of e c is not valid in those cases, but 
rather that the numerical simulations have not been per- 
formed using a transient long enough. To show the influ- 
ence of the transient time in the results, we show in Fig. 
[IJb) the dependence of e num with a for a fixed lattice of 
N = 129 sites by changing the parameter w. By increas- 
ing the transient time the numerically obtained values 
for the synchronization threshold agree better with those 
derived from transversal linear stability. The same con- 
clusions were obtained using other lattice sizes as well. 
These results suggest that the analytical result for e c re- 
mains valid, as long as we use sufficiently long transient 
times, in contrast with Ref. (l3j . 

Another factor that affects the accuracy of numeri- 
cal results for the threshold of synchronization is that 
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a distribution of initial conditions over the interval [0, 1) 
should respect the natural measure of the chaotic orbit. 
While for integer values of f3 the natural measure is uni- 
form, this is no longer valid for fractional (3, and small 
errors are introduced if we choose initial conditions with 
a uniform probability distribution. In order to overcome 
this problem we iterated each map s times before start- 
ing coupling them according to Eq. (TTJ) (this transient 
time should not be confused with the transient time T w 
we compute after having started coupling the maps). In 
Figure [He) we compare the results of two simulations: 
for the line with filled triangles we used initial conditions 
uniformly distributed along [0, 1), without discarding any 
transients (s = 0); whereas the line with open triangles 
was obtained from initial conditions chosen with respect 
to a numerical approximation of the natural measure, 
the latter having being obtained from a transient time of 
s = 10 4 iterations. The results obeying the natural mea- 
sure are more likely to agree with the analytical results 
since, after the synchronized state sets in, the correspond- 
ing orbit must follow this natural measure. 

Numerical results on the synchronization of coupled 
generalized Bernoulli maps depend also on the distribu- 
tion of initial conditions on [0,1). In order to analyze 
the dependence of e num on the initial conditions, we per- 
formed extensive numerical simulations with an ensemble 
of 5000 identically prepared CML's, each of them with a 
different initial condition [Fig. [TJd)] and the same tran- 
sient time (w = 5). We observed different values of £ num 
for each initial condition, provided a is large enough. 
Instead of showing each of them (what would turn the 
figure too much loaded with symbols) we represented in 
Fig. [ljd) only those numerical values of e nurn that are 
closest (open circles) and farthest (filled circles) with re- 
spect to the analytical value (full line). 

It is possible to understand, from a general point of 
view, the causes of the strong dependence of the synchro- 
nization threshold results on the transient time and the 
initial conditions. These causes are not restricted to cou- 
pled piecewise-linear maps as ours, but are rather generic 
for hyperbolic CML's [2fJ. We can extend our conclu- 
sions to a CML given by Eq. |T]) where the coupling pre- 
scription keeps invariant the phase space f2 = [0, 1) N , and 
for which S = {x^ = • • ■ = x^ N '} is the one-dimensional 
invariant synchronization manifold defined by the corre- 
sponding state. We consider a A-neighborhood of S as 
the set of points whose distances from the S do not ex- 
ceed A: Ea = {x : d(x, <S) < A}, where d is a suit- 
ably defined distance on the metric space f2. We define 
£ = limA— >o 5] a as a linear neighborhood of S. Accord- 
ingly r = £1 — £ is the phase space region, except the 
linear neighborhood of the synchronization manifold. 

We can speak of the global dynamics generated by the 
coupled map lattice x n+ i = BF(x„) in terms of their pe- 
riodic points. In this spirit we denote x,,- (p) the jth fixed 
point of the p-times iterated vector function BF(x„). 
The ith eigenvalue of the Jacobian matrix of BF^(x„), 
evaluated at this point, is written as Ai(xj(p)), such that 



|A 1 (x J -(p))|>--->|A iV (x J -(p))|. 

Let us consider a subset of the phase space, 4CS1, 
with natural measure fi(A). Note that, by construction, 
we have /x(£2) = 1. For hyperbolic systems satisfying the 
Axiom-A [2l| the natural measure of such subset can be 
obtained from the unstable periodic points embedded in 
it as 



14] 



tx(A) = lim VLjipl 



p— >oo 



(7) 



where Lj(p) — ]~Ii=i l^i( x i(p))l (d u 1S the largest inte- 
ger such that |Ad u (x :) (p))| > 1) and the sum sweeps over 
all Xj(p) G A. The exploitation of this identity is the 
object of periodic-orbit theory, that has been used for 
a number of theoretical investigations on the properties 
of chaotic dynamical systems [lfl [JJj]. For generalized 
Bernoulli maps (3x (mod 1) and a linear coupling, the 
Jacobian matrix has constant entries and thus do not 
depend on the orbit points, i.e. all the unstable pe- 
riodic orbits have the same eigenvalue spectra (conse- 
quently Lj(p) — L(p) for all j), and the natural measure 
is fi(A) = linip^oo Na(p)/ L(p), where Na{p) is the num- 
ber of period-p points contained in the subset A of f2. 

A byproduct of the periodic-orbit theory is that the 
(linear) transversal stability of the synchronization man- 
ifold can be studied either from the natural measure of 
a typical chaotic orbit (by the second largest Lyapunov 
exponent) or from the atypical measure generated by the 
unstable periodic orbits. In particular, with respect to 
the period-p orbit the threshold of transversal stability of 
the synchronization manifold can be obtained from the 
condition |A 2 (xj(p))| = 1 for all x,,(p) e <S. As the period 
p goes to infinity we expect an increasingly better agree- 
ment of this result with that obtained by using the second 
largest transversal Lyapunov exponent (or A2 =0). For 
a given a and values of the coupling strengths such that 
£io{oi) < £ < s up (a), that the natural measure of the 
subset A is 



li(A) = lim N A (p)/(3 p . 



(8) 



Taking A to be the linear neighborhood of the syn- 
chronization manifold, S, there follows that the number 
of orbits in this neighborhood is Ns = (3 P — 1 for integer 
(3 (if P is fractional, as in the numerical simulations of 
the previous section, one has to take the integer part of 
this expression) and the corresponding natural measure 
is given by 



fx{E) = lim (/3 P - 1)/T 



1. 



(9) 



demonstrating that the linear neighborhood of the syn- 
chronization manifold S is the asymptotic state of any 
typical initial condition (in the sense that the set of initial 
conditions that do not converge to E has zero Lebesgue 
measure). An immediate consequence of this result is 
that the natural measure outside the linear neighborhood 
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is zero since, using the fact that the natural measure is 
ergodic, we have fi(T) — — = 0. 

Given that almost all initial conditions outside the syn- 
chronization manifold eventually asymptote to it, we may 
well ask why sometimes it takes so long for this conver- 
gence to be observed in numerical experiments. As we 
saw in the previous section, this long transient time may 
even be mistaken as a effect of non-convergence. The an- 
swer lies in the properties of the horseshoe-like invariant 
chaotic set embedded in T. This set is non-attracting 
since almost all initial conditions in T converge to S as 
the time goes to infinity. 

However, the measure generated by chaotic orbits 
whose initial conditions are uniformly distributed over 
of an open region B of the phase space decays expo- 
nentially with time with escape rate 7, according to 

lime" 7 *=lim V l/Lj(p), (10) 
x 3 ( P )es 

where in the sum we consider only the periodic orbits 
of the horseshoe-like set B outside the synchronization 
manifold [l4l ]. Hence, if one picks up at random an initial 
condition off the synchronization manifold, the distribu- 
tion of the transient times is likely to be exponential, 



with a characteristic exponent dependent on the escape 
rate 7. If the initial condition is too close to an invari- 
ant manifold of an unstable periodic orbit of B it would 
stick to the manifold for some time-span and hence it 
takes a very long time for such a trajectory to approach 
the synchronization manifold. This seems to occur very 
often if we use fractional values of /3, like in the numerical 
simulations we shown in this letter. 

In conclusion, the analytical conditions for the thresh- 
old of transversal stability of the synchronized state of 
coupled piecewise-linear maps are confirmed by numeri- 
cal experiments as long as we observe the following pre- 
cautions: (i) the transient time should be chosen as large 
as possible, (ii) the choice of initial conditions should be 
done using a probability distribution which best matches 
the natural measure of the distribution. 

Although these computational problems are less likely 
to occur in coupled smooth maps, they do not invalidate 
the analytical approach to the transversal stability of 
coupled non-smooth maps, like piecewise-linear ones. We 
have used general arguments valid for hyperbolic CML's 
so as to prove that the transversal stability of the syn- 
chronized state actually implies the synchronization of 
all typical orbits. 
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